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ABSTRACT 

Techniques  for  nondestructive  failure  prediction  in  ceramics  are  examined 
in  the  context  of  a probabilistic  framework  for  obtaining  failure  and 
rejection  probabilities.  The  ultrasonic  method  appears  to  have  the  greatest 
short-term  potential  for  achieving  acceptable  failure  probabilities,  without 
excessive  rejection  of  satisfactory  components. 


Rockwell  International 
Science  Center 


SC5064.3FR 


INTRODUCTION 

Several  unique  properties  of  ceramics  are  presently  being  exploited  in 
structural  applications.  Specifically,  their  good  high  temperature  strength 
and  oxidation  (corrosion)  resistance  has  engendered  considerable  interest  in 
ceramic  energy  conversion  systems  (turbines,  heat  exchangers),  while  their 
high  hardness  renders  them  uniquely  applicable  to  bearings,  valves,  etc.  One 
of  the  primary  problems  associated  with  the  structural  application  of  such 
brittle  materials  is  their  wide  fracture  strength  variability.  In  consequence, 
effective  failure  prediction  techniques  are  a prerequisite  to  the  reliable 
use  of  ceramics  in  structural  applications. 

The  prediction  of  failure  in  structural  ceramics  can,  in  principle,  be 
achieved  using  a number  of  independent  techniques.  The  most  extensively 
developed  are  the  direct  defect  detection  techniques  (x-radiography,  pene- 
trants, and  ultrasonics)  and  the  flaw  strength  characterization  techniques 
(overload  proof  testing  and  statistical  analysis).  Many  other  methods  are 
available  as  possible  (but  less  likely)  failure  prediction  candidates;  some 
of  these  are  based  on  quite  independent  physical  phenomena,  e.g.  acoustic 
emission,  exo-electron  emission,  while  others  are  alternate  methods  of 
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visualizing  similar  physical  effects,  e.g.  acoustic  holography.  All  of  the 
available  techniques  have  limitations  when  applied  to  ceramic  systems,  and  one 
method  is  unlikely  to  emerge  as  the  universal  failure  prediction  approach. 

The  pertinent  method  will  probably  depend  on  the  component  geometry,  the 
in-service  stress  state,  and  the  material  microstructure.  Also,  several 
independent  techniques  may  be  required  to  assure  the  structural  integrity  of  a 
given  component,  without  the  excessive  rejection  of  satisfactory  parts. 

In  this  review,  an  approach  for  nondestructive  failure  prediction  in 
ceramic  systems  is  developed,  to  provide  a basis  for  assessing  the 
capabilities  of  the  potential  failure  prediction  approaches.  Then,  the  merits 
and  limitations  of  the  principal  avialable  techniques  for  failure  prediction 
are  reviewed. 


ACCEPT/REJECT  CRITERIA 

The  fracture  of  ceramics  generally  occurs  by  the  direct  extension  of  small 
cracks  in  the  vicinity  of  pre-existing  defects;  these  defects  include  large 
pores,  inclusions  and  large  grains  (frequently  as  sites  for  machining  flaws). 
The  fracture  condition  is  not,  in  general,  related  uniquely  to  the  defect 
dimensions.  A probability  of  fracture  (i(c^/a,b,c)  must,  therefore,  be  assigned 
to  a defect  of  given  size  (with  three  principal  dimensions  a,b,c)  for  a 
specific  level  of  applied  stress  cr  in  the  volume  element  containing  the 
defect  , 

The  nondestructive  capability  for  determining  the  dimensions  of  a defect 
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depends  on  the  inspection  method;  invariably,  an  error  is  associated  with  the 


estimate  of  these  dimensions.  The  magnitude  of  the  error  and  its  specific 


form  are  sensitively  dependent  on  the  method  of  inspection.  In  general,  there 
will  be  a specific  probability, e. g.  <6(agj/a),  for  each  technique  that  the 
estimated  defect  dimensions  b^^,  Cg^jwill  be  in  a certain  size 

range,  given  the  actual  defect  dimensions  (a,b,c). 


We  have  thus  far  identified  two  probabilities:  one  concering  the  fracture 


probability  for  a defect  of  given  dimensions,  the  other  concering  the 


nondestructive  estimate  of  the  defect  dimensions.  However,  a third 


probability  is  also  needed  to  complete  the  reliability  analysis;  namely,  the 


probability  that  a defect  in  a given  size  range  will  exist  in  the  volume  of 


the  component  being  inspected,  tCa).  This  probability  is  estimated  by 


metal lographic  studies  on  samples  taken  from  each  batch  of  material. 


The  product  of  these  three  probabilities  can  be  integrated  to  various 


inspection  levels,  a^^,  to  obtain  two  interrelated  probabilities:  the 


false-accept  probability and  the  false-reject  probability  (Fig.  la) 


*A  ' jT  J f [<i>(of/a)<lo-,jj<f)(ag5/a)dagJ[<i>(a)daJ 


0 0 0 
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STRESS  LEVEL,  a. 


INSPECTION  LEVEL,  b? 


fail  in  service,  i.e.  is  the  in-service  failure  probability.  This 
probability  decreases,  of  course,  as  the  inspection  level  decreases 
(Fig.  la).  The  false-reject  probability  is  the  (related)  probability 
that  rejected  components  would,  in  fact,  have  performed  sat’ sfm  tori ly  in 

service,  because  the  defect  severity  has  been  overestimated  by  the  selected 

★ 

inspection  level.  This  probility  increases  as  a^^  decreases  (Fig.  la). 
However,  it  is  crucial  to  recognize  that  these  probabilities  are  interrelated, 
i.e.  they  merely  represent  different  ranges  of  integration  of  the  same 
combination  of  probability  functions  (Eq.  1).  This  interdependence  is 
exemplified  in  Fig.  lb,  which  is  a typical  plot  relating  the  false-accept  and 
false-reject  probabilities  - once  one  of  these  probabilities  has  been 
selected,  the  other  probability  as  well  as  the  associated  inspection  level  are 
necessarily  defined.  It  is  now  apparent  from  Fig.  lb  that  the  inspection 
technique,  or  combination  of  techniques,  that  would  be  preferred  is  that  which 
yields  a curve  as  close  as  possible  to  the  probability  axes.  For  example,  in 
Fig.  lb,  technique  B is  preferred  over  technique  A,  because  the  rejection  cf 
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satisfactory  components  required  to  satisfy  the  failure  probability 
requirement  is  much  lower.  Such  curves  thus  represent  a quantitative  method 
for  characterizing  the  failure  prediction  capabilities  of  various  inspection 
techniques,  for  a given  material  and  service  condition.  However,  the 
generation  of  these  curves  represents  a major  experimental  effort,  and  the 
information  presently  available  is  limited.  The  present  status  and  imminent 
developments  are  examined  in  subsequent  sections.  It  is  hoped  that  all  future 
studies  will  be  directed  toward  the  development  of  false-accept,  false-reject 
curves,  to  enable  a quantitative  framework  for  failure  prediction  in 
structural  ceramics  to  be  constructed. 

failure  models 

The  fracture  process  in  a ceramic  depends  specifically  on  the  material  and 
the  defect  responsible  for  fracture.  In  most  materials,  surface  cracks 
introduced  by  machining  are  a consistent  source  of  failure.  This  failure 
mechanism  is  particularly  prevalent  in  dense,  coarse  grained  materials. 
Inclusions  introduced  during  fabrication  are  another  prominent  origin  of 
failure,  particularly  in  dense,  fine  grained  materials  (for  which  failure  from 
surface  cracks  is  less  dominant).  Individual  large  pores,  or  a concentrated 
array  of  fine  pores,  are  another  frequently  observed  source  of  failure, 
particularly  in  sintered  materials.  Finally,  surface  cracks  introduced  by 
post  fabrication,  environmental  interaction  phenomena  can  be  important,  i.e. 
cracks  introduced  by  projectile  impact,  contact  stresses  etc.  or  cracks  formed 
by  oxidation. 
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Probablistic  infonnation  concerning  these  failure  processes  is  limited. 

The  most  extensive  information  is  available  for  hot  pressed  silicon  nitride. 

This  material  is  thus  used  as  an  example  of  the  probablistic  treatment  of 
failure.  However,  it  is  re-emphasized  that  the  failure  model  can  be  highly 
specific,  both  to  the  host  material  and  the  defect  type.  Therefore,  the 
models  described  should  only  be  extended  to  other  materials/defects  with 
considerable  caution. 

Surface  Cracks 

The  surface  crack  has  been  the  subject  of  extensive  recent  analysis.  In 
consequence,  the  stress  intensity  factor  distribution  is  now  relatively  well 
comprehended  for  semi-circular  cracks  normal  to  the  surface  (Fig.  2)  and 
reasonably  we  11 -character! zed  for  normal  semi -elliptical  cracks.  The 
equivalent  solutions  for  cracks  at  other  inclinations  to  the  surface  are  not 
known  in  detail;  although  reasonable  estimates  can  be  based  on  results  for  | 

contained  elliptical  cracks.  | 

One  possible  model  of  fracture  from  surface  cracks  considers  that  fracture 

I 

occurs  when  the  peak  value  of  the  stress  intensity  factor  reaches  the  local 
toughness  of  the  material,  yielding; 

(7-^  * Z (a/c,fl)  Kj.  a"^  (2) 

where  cr^  is  the  fracture  stress  measured  in  the  absence  of  slow  crack  I 

growth,  a and  c are  the  two  principal  dimensions  of  the  crack,  9 is  its 
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inclination  to  the  surface,  and  Z is  the  analytically  determined  stress 
intensity  function.  For  this  model,  the  variability  in  the  fracture  stress  is 
attributed  exclusively  to  variations  in  the  local  toughness:  distributions  in 
oj,  (<i((r^/a,c,d))  pertinent  to  such  variability  are  likely  to  be  Gaussian. 
However,  sub-critical  extension  of  the  crack  in  a region  of  high  local  stress 
intensity  factor  (Fig.  2)  may  precede  final  fracture.  Then,  there  may  be 
systematic  trends  in  the  fracture  stress  with  a/c  and  d . Probability 
functions  that  neglect  these  possible  trends  will  exhibit  a wider  and  less 
characteristic  distribution  than  those  that  take  these  effects  into  account. 
The  details  of  the  fracture  model  thus  constitutes  an  important  part  of  the 
reliability  analysis. 

Specific  probability  functions  for  surface  cracks  have  not  yet  been 
obtained,  largely  because  of  the  difficulty  in  distinguishing  the  profile  of 
the  fracture  initiating  crack  on  the  fracture  surface.  A viable  solution  to 
this  difficulty  has  not  yet  been  identified.  Some  useful  information  can  be 
obtained  by  introducing  well-defined  cracks  using  the  Knoop  indentation 
technique  (Fig.  3).  However,  this  method  samples  only  the  variability  in  the 
fracture  toughness.  Much  additional  work  is  needed  in  this  area. 

Voids 

Fracture  from  isolated  voids  in  ceramics  has  been  analyzed  using  several 
models.  The  model  which  accounts  most  effectively  for  the  observed  trends 
with  void  size  is  a statistical  model,  which  considers  that  fracture  occurs  by 
the  activation  of  a distribution  of  flaws  that  pre-exist  in  the  vicinity  of 
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the  void  (Fig.  1).  The  most  likely  fracture  initiating  flaws  are  microcracks 
associated  with  grain  boundary  cusps  on  the  void  surface  (small  inclusions 
near  the  void  are  an  alternate  possibility).  Analyses  of  both  surface  and 
volume  distributed  flaws  have  been  conducted  by  considering  a flaw  strength 
distribution  of  the  type 

g(S)  «A(S/Sq)''  (3) 

where  g(S)dS  is  the  number  of  flaws  in  unit  area  (volume)  with  a strength 
(i.e.  extension  stress  in  uniaxial  tension)  between  S and  S+dS,  and 
and  k are  distribution  parameters  to  be  determined  from  the  test  data.  For 
the  surface  flaw  problem,  the  analysis  encompasses  all  possible  flaw  sizes, 
and  yields  the  general  result  that  the  probability  of  fracture  ♦(o’f)  at  an 
applied  stress  is  given  by; 

D*'(a)  exp  (0.52k  - 0.14)  (4) 

2 

where  r is  the  void  radius,  a*  (1/r) (K^/S^^)  and  D(ck)  is  the  function 
plotted  in  Fig.  4.  For  the  volume  flaw  problem,  the  analysis  has  only  been 
conducted  for  flaws  very  much  smaller  than  the  void  radius,  and  is  thus 
restricted  in  utility.  Additional  work  is  needed  to  extend  the  analysis  into 
the  large  flaw  regime. 

A considerable  quantity  of  fracture  data  are  required  to  determine  the 
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detailed  applicability  of  the  surface  model  vis-a-vis  the  volume  model  for  a 
particular  material.  Sufficient  data  for  this  purpose  are  not  presently 
available.  However,  the  general  pertinence  of  a statistically  based  model 
seems  to  be  substantiated  by  the  available  data,  particularly  by  the  size 
dependence  of  the  fracture  stress  at  constant  probability:  other  models 
predict  either  no  size  dependence  (stress  concentration  model)  or  on  inverse 
square  root  size  dependence  (crack  equivalence  models).  Values  of  the 
distribution  parameters  can  be  obtained  from  the  data. 

Inclusions 

The  pertinent  inclusion  fracture  model  depends  on  the  specific  inclusion. 
Three  types  of  inclusions  in  silicon  nitride  are  considered,  to  illustrate  the 
range  of  possibilities:  silicon,  tungsten  carbide,  and  iron  silicide.  The 
trends  with  inclusion  size  obtained  for  these  inclusions  are  depicted  in 
Fig.  5.  The  silicon  inclusion  has  a strong  influence  on  the  strength;  whereas 
the  effect  of  the  tungsten  carbide  inclusion  is  minimal.  The  need  to 
nondestructive ly  identify  the  defect  type  is  thus  dramatically  demonstrated. 

Post-fracture  inspection  of  the  silicon  inclusions  indicates  that  they  are 
dense., i.e.  contain  little  porosity,  and  exhibit  transgranular  fracture.  These 
observations,  coupled  with  the  low  toughness  of  silicon  (0.6  MPavTii),  suggest 
that  the  silicon  inclusions  fracture  sub-critically,  i.e.  without  propagating 
into  the  much  tougher  ( 5 MPa  vIS)  silicon  nitride  matrix.  This  expectation  is 
verified  by  acoustic  emission  studies  on  stressed  samples  containing  silicon 
inclusions.  Hence,  since  the  elastic  properties  of  silicon  are  similar  to 
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those  of  silicon  nitride,  a plausible  model  of  ultimate  fracture  considers  the 
inclusions  as  cracks  of  equivalent  dimensions  {Fig.  6).  Then,  the  fracture 
stress  is  governed  by  an  equation  similar  to  Eq.(2)  that  contains  the 
appropriate  Z(a/c)  for  an  internal  crack.  A comparison  of  measured  strengths 
with  strengths  predicted  by  the  crack  equivalence  model  for  an  assumed 
toughness  of  5 MPa-Zm  is  presented  in  Fig.  6.  The  variability  in  the 
measured  fracture  stress  might  be  attributed  either  to  a toughness  variability 
or  to  a combined  toughness,  sub-critical  extension  effect.  The  data  are  not 
comprehensive  enough  to  unequivocally  separate  these  possibilities,  but  it  is 
noted  that  the  ratio  of  the  measured  to  the  predicted  strength  exhibits  a 
close  conformance  to  a normal  distribution.  This  type  of  distribution  would 
be  anticipated  for  the  toughness  when,  as  in  the  present  case,  the  crack  front 
samples  a large  number  of  silicon  nitride  grains.  The  center  of  this 
distribution  at  0.8  indicates  either  that  the  average  local  toughness  is  less 
then  the  assumed  value  (4  compared  with  5 MPa-Zm),  or  that  the  crack 
dimensions  at  criticality  are  consistently  larger  than  the  inclusion 
dimensions,  because  of  sub-critical  extension  phenomena.  Finally,  note  that 
the  important  inclusion  dimensions  are  the  two  dimensions  (a,c)  normal  to  the 
applied  tension. 

Inspection  of  the  fractured  iron  silicide  inclusions  indicates  extensive 
porosity  and  an  irregular  fracture  surface.  A different  fracture  mechanism 
might  thus  pertain.  This  is  verified  when  the  above  fracture  model  is  applied 
to  test  data;  the  measured  strengths  greatly  exceed  the  predicted  values  and 
exhibit  a diminished  dependence  on  the  inclusion  size.  An  alternative  model 
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suggested  by  the  fracture  morphology  is  that  fracture  initiates  from  the  pores 
within  the  inclusions,  i.e.  the  inclusion  is  treated  as  a porous  body 
subjected  to  the  stress  field  imposed  through  the  silicon  nitride  matrix. 

This  stress  field  is  the  result  of  thermal  expansion  mismatch  and  the  applied 
stress.  For  ellipsoidal  inclusions,  both  of  these  stresses  are  hydrostatic 
(equal  in  all  three  directions),  and  the  analysis  of  fracture  then  becomes  a 
simple  statistical  problem.  Analyzing  the  porous  inclusion  using  statistical 
flaw  distributions,  in  the  manner  described  above  for  fracture  from  voids,  the 
probability  of  fracture  becomes 


♦ * 1-exp 


0^  -t-q-g 


k 


(5) 


where  V is  the  volume  of  the  inclusion,  is  the  thermal  mismatch  stress  and 
/3  is  the  ratio  of  the  applied  stress  to  the  stress  in  the  inclusion.  Test 
data  analyzed  according  to  this  model  are  plotted  in  Fig.  7.  Reasonable 
values  of  the  parameters  are  obtained,  indicating  that  the  model  has  merit. 
Note,  therefore,  that  the  inclusion  parameter  of  interest  in  this  instance  is 
its  total  volume,  V. 

The  minor  effect  of  tungsten  carbide  inclusions  on  the  strength  derives 
from  both  their  high  toughness  and  large  elastic  modulus.  The  high  toughness 
prevents  premature  inclusion  fracture,  while  the  large  modulus  confines  the 
zone  of  tension  in  the  matrix  to  a small  region  around  the  poles  of  the 
inclusion.  In  consequence,  the  statistical  analysis  of  fracture  from  this 
small  zone  of  tension  indicates  low  fracture  probabilities  at  appreciable 
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levels  of  the  applied  stress.  In  qualitative  confirmation  of  these 
expectations,  the  fracture  of  samples  containing  tungsten  carbide  inclusions 
usually  occurs  from  small  matrix  voids  in  the  vicinity  of  the  poles.  There 
are  not  yet  sufficient  fracture  data  to  identify  the  distribution  parameters 
pertinent  to  this  fracture  process;  but,  tentatively,  the  fracture  probability 
for  these  inclusions  can  be  equated  to  zero,  for  the  stress  levels  of 
practical  interest,  ^ 300  MPa  (Fig.  5). 

DEFECT  SIZE  ESTIMATION 

It  Is  apparent  from  the  fracture  models  that  greatly  reduced  false-reject 
probabilities  will  result  if  the  defect  type,  as  well  as  its  pertinent 
difflension(s),  can  be  elucidated  by  the  inspection  method.  This  requirement  is 
a vital  consideration  in  the  comparison  of  the  various  methods  of  defect 
determination.  This  particular  problem  is  not  a concern  for  the  indirect 
failure  prediction  method,  which  estimates  the  defect  'strength';  but  other 
problems  are  introduced.  The  indirect  methods  are  thus  considered  first, 
followed  by  a discussion  of  the  direct  defect  detection  methods. 

Indirect  Methods 

Statistical  Analysis 

The  statistical  prediction  of  failure  relies  on  an  ability  to  directly 
characterize  the  flaw  strength  distribution  function.  The  strength 
distribution  function  g(S)dS  is  defined  as  the  number  of  flaws  per  unit  volume 
(or  unit  area  for  surface  flaws)  with  a strength  between  S and  S+dS.  For 
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non-interacting  flaws,  this  function  is  related  to  the  fracture  probability 
#(S)  by 


♦(S)  » 1 -exp 


(6) 


where  is  the  sample  volume.  Recently,  a method  for  obtaining  g(S)  from 
fracture  data  has  been  devised.  A typical  example,  which  applies  to  tensile 
tests  is 


g(S)  « t/Trrh  (7) 

where  r is  the  specimen  radius,  1 is  the  gauge  length,  and  -ln(l-*(S)). 
Once  g(S)  has  been  obtained  in  this  way,  within  the  stress  range  that  exists 
in  the  test  component,  Eq.(6)  can  be  evaluated  numerically  for  any  specimen 
geometry  and  stress  distribution.  The  constraint  that  the  test  data  pertinent 
to  low  failure  probabilities  be  obtained  in  the  stress  range  of  concern 
(because  extrapolation  to  low  stresses  cannot  be  performed  with  appreciable 
confidence)  means,  for  most  applications,  that  extensive  tensile  data  on 
relatively  large  samples  are  required.  (By  contrast,  the  low  failure 
probabilities  for  the  direct  methods  relate  primarily  to  the  probability, 
<6(a)da  (see  Eqn.  1),  which  can  be  determined  metallographically.) 

The  method  in  its  present  form  has  an  additional  weakness  - notably,  the 
requirement  that  the  flaw  population  be  invariant  from  component  to 
component.  Such  an  invariance  is  quite  unlikely  for  typical  ceramic 
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fabrication  processes  and  hence,  high  confidence  in  the  predicted  fracture 
probability  cannot  be  achieved.  Some  independent  method  for  identifying 
variations  in  the  flaw  population  is  thus  required  before  statistical  methods 
can  be  realistically  applied.  Such  a possibility  has  recently  been  proposed, 
based  on  the  ultrasonic  attenuation  a of  the  component. 

Ultrasonic  attenuation  affords  a measure  of  the  large  extreme  of  the 
microstructure,  or  of  the  surface  crack  population; 


a(f) 


/«(  r,f)g(r)dr 
0 


(8) 


where  f is  the  frequency,  Q is  the  cross  section  of  the  scatterer,  and  g(r)dr 
is  the  number  of  large  grains  (pores,  inclusions,  surface  cracks)  per  unit 
volume  (area)  in  the  size  range  r to  r+dr.  Hence  for  materials  in  which 
fracture  is  determined  by  the  large  extreme  of  the  microstructure  (e.g.  large 
pores  in  reaction  bonded  silicon  nitride  and  large  grains  in  relatively  coarse 
grained  materials),  or  by  surface  cracks,  g(r)dr  is  related  to  the  flaw 
strength  distribution  g(S)dS  through  the  usual  relation 


''c 

TTT 


(9) 


where  Y is  a parameter  that  depends  on  the  flaw  shape.  Hence,  an  attenuation 
measurement  could  define  g(S)dS  and  thereby  the  fracture  probability  #(S). 
Variability  in  the  flaw  population  from  component  to  component  might  thus  be 
directly  incorporated  into  the  failure  prediction.  This  intriguing 
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Overload  Proof  Testing 

The  theory  of  proof  testing  has  been  established  for  several  years,  but 

the  pertinent  experimental  evaluation  is  only  now  in  progress.  The  present 

status  is  as  follows.  The  utility  of  overload  proof  testing  to  eliminate 

those  components  that  would  normally  fail  during  service  has  been  established 

theoretically  and  experimentally  for  proof  tests  conducted  in  an  ambience  that 

excludes  significant  slow  crack  growth  during  the  test.  For  such  conditions 

the  minimum  time  to  failure,  t.^.,  is 

wm 


min 


(10) 


where  R is  the  proof  ratio  and  n,  v^^  and  are  slow  crack  growth 
parameters  (v  » There  is  also  a well-defined  variance  of 

the  failure  time  that  depends  on  the  variances  in  the  slow  crack  growth 
parameters  and  the  proof  ratio.  This  minimum  failure  time  is  realized 
whenever  the  stress  in  each  element  of  the  component  during  the  proof  test 
exceeds  the  service  stress  (in  the  equivalent  element)  by  an  amount  greater 
than  R;  provided  that  extraneous  interface  stresses  are  not  developed  during 
unloading  and  that  additional  defects  are  not  produced  (e.g.  by  projectile 
impact  and  oxidation)  during  service. 

When  it  is  not  possible  to  exclude  slow  crack  growth  during  the  proof 
test,  theory  indicates  that  a minimum  unloading  time  is  needed  to  assure  a 
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finite  survival  time  after  proof  testing.  However,  even  if  this  minimum  is 
impractical,  the  theory  predicts  that  the  failure  probability  after  proof 
testing  can  be  many  orders  of  magnitude  less  than  the  failure  probability 
without  proof  testing.  This  latter  effect  has  been  observed,  but  the  failure 
probability  is  often  larger  than  the  predicted  value.  This  disparity  is 
probably  associated  with  the  presence  of  crack  growth  instabilities  that  occur 
during  unloading  at  stress  intensities  close  to  the  critical  value. 

Another  restriction  of  proof  testing  concerns  the  practicality  of 

reproducing  the  in-service  stress  distribution  in  the  proof  test.  Except  for 

simple  geometries  or  service  environments,  a substantial  fraction  of  volume 
elements  in  the  proof  test  will  experience  stress  ratios  that  differ  from  the 

ideal  value  specified  by  Eq.{10).  Those  elements  for  which  R is  smaller  than 

ideal  introduce  a false-accept  probability,  and  those  for  which  R is  greater 
than  ideal  introduces  a false-reject  probability.  Hence,  each  proof  test  can 
be  expressed  as  a false-accept,  false-reject  curve  for  comparison  with  the 
curves  derived  for  the  direct  methods.  This  quantification  of  proof  tests  has 
not  yet  been  implemented,  but  it  is  urged  that  this  method  of  gauging  the 
relative  efficacy  of  a particular  proof  test  design,  vis-a-vis  the  direct 
nondestructive  methods,  be  adopted  as  a future  requirement. 

Direct  Methods 

A wide  range  of  methods  for  direct  defect  detection  has  been  devised 
(Table  I).  However,  only  the  ultrasonic  and  x-ray  methods  have  been 
sufficiently  developed  that  some  of  the  important  issues  concerning  defect 
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TABLE  I 


INSPECTION  METHODS  PRIMARY  SENSITIVITY 

BULK  DEFECTS 


1. 

X-Radlography  (Microfocus) 

Density  (Atomic  Number) 

2. 

Ultrasonics 

Acoustic  Impedance 

NEAR  SURFACE  DEFECTS 

1. 

Ultrasonics  (Surface  Waves) 

Acoustic  Impedence 

2. 

Microwaves 

Dielectric  Constant 

3. 

Acoustic  Microscopy 

Acoustic  Impedance 

NOVEL  APPROACHES 

Photo-Acoustic  Spectroscopy  (Optical  Absorption) 
Eddy  Currents  (Space  Charge  Interaction) 
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distinction  and  size  characterization  can  be  adequately  discussed.  The 
detailed  dicussion  is  restricted  to  these  two  methods,  but  this  should  not  be 
construed  to  indicate  that  the  other  methods,  as  they  develop,  will  not 
ultimately  play  an  important  role  in  the  overall  failure  prediction  scheme. 

For  example,  the  ability  of  microwave  methods  to  detect  silicon  inclusions 
might  be  of  great  significance  in  failure  prediction  in  silicon  nitride;  also, 
improved  penetrant  techniques  (e.g.  dye  enhanced  x-ray  methods)  may  perhaps 
find  an  important  niche  in  surface  crack  characterization  schemes. 


X-Radioqraphy 

X-radiography  is  a we 11 -developed  technique,  and  the  principles  involved 
in  the  generation  of  images  from  defects  are  largely  understood.  A brief 
statement  of  these  principles  allows  judgments  to  be  made  about  the  role  of 
X-radiography  in  ceramic  failure  prediction.  The  two  image  variables  are  the 
resolution  and  the  image  contrast.  The  resolution  x is  related  to  the  focal 
spot  diameter  f,  and  the  distances  d^,  between  the  focal  point  and  the 
defect  and  d2,  between  the  defect  and  film  by 

Conventional  x-ray  equipment  utilizes  focal  spot  tubes  ^SOO^m  in  diameter  and 
hence,  for  typical  dispositions  of  the  source,  component  and  film  (as  realized 
for  nonsimple  component  geometries)  the  maximum  resolution  (:5100»xm)  is 
inadequate  for  most  structural  ceramic  applications.  However,  refined  x-ray 
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systems  are  capable  of  achieving  much  smaller  focal  spot  sizes;  in  fact,  the 
Cosslett-Nixon  type  projection  instrument  (which  utilizes  electron  focussing 
tubes)  is  capable  of  producing  spot  sizes  as  small  as  ^iMm.  Hence,  resolution 
is  not  a limitation  per  se  to  the  application  of  X-radiography  to  structural 
ceramics. 

The  image  contrast  represents  a more  serious  limitation.  The  contrast, 

AI,  on  the  film  is  given  by 

Al  « 0.86am  aV  q//3  (12) 

where  Am  is  the  difference  in  linear  absorption  coefficient  between  the  defect 
and  the  matrix,  2a  is  the  defect  thickness  (in  the  direction  of  the  beam), 

Yjj  is  the  x-ray  "ganma"  of  the  film  and  0 is  a numerical  build-up  factor. 

The  linear  absorption  coefficient  is  primarily  determined  by  the  atomic  number 
and  density  of  the  material  and  varies  with  the  energy  of  the  incident  x-ray 
photons.  The  atomic  number/density  relation  for  the  elements  with  closely 
similar  atomic  numbers  can  exhibit  large  x-ray  absorption  differences  and  vice 
versa.  Certain  types  of  defects  will  thus  be  capable  of  generating 
significant  image  contrast  for  a given  matrix  material;  these  defect  types  can 
often  be  predicted  from  available  absorption  tables.  The  absorption 
differential, Am,  for  a given  defect  Increases  quite  rapidly  as  the  x-ray 
photon  energy  decreases  and  hence,  to  maximize  detectability,  the  minimum 
energy  consistent  with  a practicable  exposure  time  should  be  selected.  The 
build-up  factor,  0,  increases  as  the  sample  thickness  increases  and  thus,  the 
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component  thickness  is  a variable  effecting  defect  detectability.  Finally, 
since  is  directly  proportional  to  the  defect  thickness  a,  there  will  exist 
a lower  limit  to  the  defect  size  that  yields  detectable  image  contrast. 

The  quantities  Tq  and  $ cannot  be  effectively  predetermined;  hence, 
empirical  guidelines  on  defect  detectability  must  firstly  be  obtained  for 
typical  defect/matrix  systems,  and  then  approximate  scaling  effects  can  be 
predicted  using  Eq.(12).  Preliminary  studies  on  inclusions  in  silicon  nitride 
indicate  that  inclusions  containing  high  atomic  number  elements  (e.g.  W,  Fe) 
are  readily  detected  using  microfocus  systems,  at  the  smallest  sizes  of 
interest  (~25Mm).  The  detectability  of  inclusions  with  comparable  atomic 
numbers  to  the  matrix,  such  as  Si  and  SiC,  voids  and  cracks  is  less 
satisfactory;  although,  the  detection  limits  that  can  be  achieved  with  image 
enhancement  capabilities  has  yet  to  be  ascertained.  Similarly,  the  ability  to 
distinguish  different  defect  types  has  yet  to  be  evaluated.  Some  information 
is  clearly  available  in  the  contrast  (Eq.  12),  but  the  extent  to  which  this 
provides  a satisfactory  classification  remains  to  be  determined.  Another 
difficulty  that  must  be  surmounted  concerns  the  ability  to  obtain  the  defect 
dimension  normal  to  the  surface  . This  dimension  Influences  the  contrast 
(Eqn.  12),  but  it  is  not  likely  to  be  separable  from  the  effect  of  the  defect 
type.  Tomographic  techniques  are  a more  plausible  probability;  but  again, 
their  efficacy  in  components  of  complex  shape  has  yet  to  be  determined. 


The  nature  of  the  probability  function  <<>(agj| a)  pertinent  to  the  x-ray 
method  deserves  brief  mention.  Since  the  method  creates  an  image  of  the 
defect,  the  error  is  concerned  both  with  the  measurement  of  the  dimensions 
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from  the  image,  and  the  ability  to  obtain  a three-dimensional  image.  The 
former  represents  a random  error,  and  is  probably  Gaussian;  while  the  latter 
is  a systematic  error  related  to  the  component  geometry  and  its  effect  on 
accessability  of  the  x-ray  beam  to  the  defect.  The  magnitudes  of  these  errors 
are  presently  being  evaluated. 

Finally,  it  is  instructive  to  cursorily  explore  the  utility  of  the  x-ray 
method  for  failure  prediction  in  hot  pressed  silicon  nitride.  If  Si 
inclusions^ SOOwm  cannot  be  discerned,  (as  the  available  information  suggests) 
then,  since  Si  inclusions  ^150/im  in  diameter  are  potential  critical  defects 
(Fig.  5),  it  may  be  concluded  that  the  x-ray  method  is  incapable  of  estimating 
the  dimensions  of  certain  types  of  critical  defects.  The  false-accept  (or 
failure)  probability  for  the  x-ray  method  is  thus  identical  to  the  probability 
of  occurrence  of  Si  inclusions,  in  the  size  range ~ 150  to  SOOMm.  The  method 
has  not,  therefore,  effected  a substantial  reduction  in  the  false-accept 
probability,  vis-a-vis  its  value  without  inspection.  This  result  suggests 
that  future  x-ray  development  pertinent  to  silicon  nitride  should  be  devoted 
to  enhancing  the  detectability  of  inclusions  of  this  type. 

Ultrasonics 

The  advances  in  ultrasonic  defect  characterization  have  been  appreciable 
in  recent  years.  It  is  now  apparent  that  all  of  the  bulk  defects  and  most  of 
the  surface  defects  of  concern  in  ceramics  can  be  detected  with  the 
appropriate  ultrasonic  technique.  However,  much  additional  study  is  needed  to 
determine  how  well  the  pertinent  defect  dimensions  and  the  defect  type  can  be 
characterized.  The  requisite  defect  detectability  has  been  achieved  by  noting 
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that  the  amplitude  of  the  ultrasonic  wave  scattered  by  the  defect  Increases 
rapidly  as  the  wavelength  \ decreases,  to~2TaA.  Detectability  is  thus 
inhanced  by  selecting  frequencies  (up  to  200  MHz)  in  excess  of  those 
conventionally  used  for  defect  detection  (up  to^20  MHz).  It  also  appears 
that  the  acoustic  impedance  mismatch  between  the  defect  and  the  matrix  is 
sufficient  for  most  defects  of  interest  (especially  when  the  defect  contains 
porosity  wnich  reduces  its  effective  impedence)  to  permit  the  required 
detectability. 

The  background  scattering  (which  often  imposes  a lower  limit  on  defect 
detectabiity)  is  not  a concern  for  bulk  defects  of  critical  size  ('•25Mm)  in 
fine-grained  structural  ceramics  (hot  pressed  pressed  and 

sintered  SiC);  because  the  grain  scattering  is  minimal  up  to '-400  MHz.  The 
background  is  a concern,  however,  for  surface  crack  detectability.  The 
background  in  this  case  is  associated  with  the  scattering  from  arrays  of 
surface  cracks  introduced  by  surface  grinding.  Typical  detectability  limits 
are  indicated  in  Fig.  8.  On  polished  surfaces,  the  smallest  cracks  that  can 
be  introduced  {60nm  diameter)  are  easily  detected  (Fig.  8a);  for  surfaces 
ground  with  304m  diamond  particles,  the  smallest  cracks  that  can  be  detected 
are~80Mm  in  diameter;  while,  for  rough  ground  surfaces,  cracks  smaller  than 
-1204ffl  cannot  be  detected. 

The  specific  ultrasonic  methods  for  bulk  and  surface  defects  are 
different,  but  the  principles  are  the  same.  Sulk  defects  are  characterized 
using  longitudinal  and/or  shear  waves  propagated  directly  into  the  component; 
preferably  using  a transducer  array  that  permits  rapid  electronic  scanning  of 
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components,  through  a ceramic  buffer  contoured  to  match  the  component. 
Alternatively,  a traversable  transducer  focused  through  a liquid  metal  medium 
may  be  employed.  Surface  cracks  are  characterized  using  surface  waves,  which 
can  be  excited  on  the  component  using  ceramic  waveguides  and  a fluid  or 
polymeric  couplant.  Additional  studies  are  needed,  however,  to  design  optimum 
transducer  configurations. 

The  ultimate  approach  for  ultrasonic  defect  characterization  in  ceramics 
remains  to  be  developed.  However,  present  knowledge  suggests  the  following 
sequence.  Following  an  initial  rapid  scan  to  determine  the  location  of 
defects,  the  defect  type  will  be  identified  from  the  characteristics  of  the 
scattered  signal  at  high  frequencies  (100-400  MHz).  The  signal  in  the  time 
domain  has  very  characteristic  features,  as  indicated  by  comparing  the  void 
signal  (Fig.  9a)  with  the  signal  from  a WC  inclusion  (Fig.  9b).  Additionally, 
the  frequency  domain  (obtained  using  a fast  Fourier  transform)  contains 
information  characteristic  of  the  defect  type.  Note,  however,  that  the 
pertinent  information  can  only  be  obtained  experimentally  if  transducer 
response  i?  subtracted,  e.g.  using  an  inverse  filter.  Finally,  the  frequency 
dependence  of  the  scattered  amplitude  in  the  long  wavelength  limit  (X?10a), 
i.e.  at  between  approximately  5 and  50  MHz,  can  be  used  to  provide  a unique 
measure  of  the  volume  of  the  defect  (Fig.  10).  The  scattering  model  used  to 
obtain  this  information  has  been  we  11 -substantiated  experimentally,  at  least 
for  defects  of  regular  shape,  (vis.  generalized  ellipsoids,  surface  cracks). 
The  model  indicates  that  the  scattered  amplitude  A at  low  frequencies  is  given 
by; 


31 

J/0791A/S.N 


Rockwell  International 

Science  Center 


SC5064.3FR 


l/a  FOR  VOID  IN  SLN4. 


Rockwell  International 

Scwnce  Center 


SC5064.3FR 


SC78-1484 


Figure  10 


34 


Rockwell  International 

Science  Center 


A « Vt?  ^£0^ 


SC5064.3FR 

(13) 


where  at  is  the  frequency,  and  17  Is  a parameter  that  depends  on  the  defect  type 
and  shape:  for  example,  17  for  a spherical  void  is. 


V * 


1 + V 

T-rn; 


+ 


10(1-21/) 


(14) 


where  u is  Poisson's  ratio  and  C|  is  the  longitudinal  wave  velocity  in  the 
matrix.  Individual  defect  dimensions  (a,b,c)  and  orientations  can  also  be 
obtained  from  long  wavelength  data,  if  information  is  available  for  the  same 
defect  as  several  (up  to  5)  scattering  angles.  The  viability  of  obtaining  the 
quantity  of  information  required  to  determine  individual  dimensions  has  yet  to 
be  demonstrated  for  a practical  system. 

The  error  associated  with  the  estimate  of  the  defect  volume  required  for 
the  reliability  analysis  has  not  been  evaluated  in  detail  for  real  defects; 
although  the  appropriate  experiments  are  now  underway.  The  errors  are 
expected  to  be:  a random  error  in  the  amplitude  associated  with  the  level  of 
the  background  (probably  Gaussian),  a systematic  error  if  the  defect  type  is 
incorrectly  assigned  from  the  high  frequency  analysis,  and  an  error  (probably 
random  in  most  cases)  associated  with  the  unknown  orientation  of  the  axes  of 
the  defect.  These  errors  combine  to  yield  the  required  function,  <6 (V^^l  V) . 


SUMMARY 

The  general  approach  for  quantitative  failure  prediction  in  ceramics  using 
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nondestructive  methods  of  defect  characterization  has  been  described.  The 
current  state  of  knowledge,  as  it  impinges  on  the  quantitative  failure 
prediction  issue,  has  been  reviewed.  The  knowledge  is  incomplete,  but 
considerable  progress  has  been  realized  in  the  last  two  years  toward  attaining 
the  requisite  information.  The  next  two  years  should  see  progress  at  an  equal 
or  more  rapid  rate,  and  the  emergence  of  the  first  quantitative  nondestructive 
failure  prediction  scheme  for  ceramics.  At  that  juncture,  it  will  be  possible 
to  select  the  best  possible  technique,  or  combination  of  techniques,  pertinent 
to  a specific  material  and  component  geometry.  Substantial  scope  for  the 
improvement  of  existing  inspection  techniques,  and  for  the  introduction  of  new 
techniques,  will  exist  at  that  stage;  such  advances  are  to  be  encouraged, 
using  the  quantitative  basis  for  technique  comparison  developed  in  this  review. 
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